Three and Four-Body Interactions in Spin-Based Quantum Computers 
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In the effort to design and to construct a quantum computer, several leading proposals make use 
of spin-based qubits. These designs generally assume that spins undergo pairwise interactions. We 
point out that, when several spins are engaged mutually in pairwise interactions, the quantitative 
strengths of the interactions can change and qualitatively new terms can arise in the Hamiltonian, 
including four-body interactions. In parameter regimes of experimental interest, these coherent 
effects are large enough to interfere with computation, and may require new error correction or 
avoidance techniques. 
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Researchers in the field of quantum computation 
seek to exploit quantum coherence to speed calculation. 
While theoretical possibilities have been enticing, the 
prospect of realizing a practical quantum computer is 
quite daunting. A number of physical systems have been 
suggested as candidate instantiations. One particularly 
promising proposed design 1] involves electron spins lo- 
calized in quantum dots. In its original form, the pro- 
posal supplements tunable exchange interactions between 
electron spins with single-qubit operations to achieve a 
universal quantum computer. Subsequent research has 
shown that a universal set of gates can be constructed 
using the exchange interaction alone, provided that one 
encodes a logical qubit into the state of several spins 
("encoded universality") 0]. Motivated by the scheme 
of Ref. l], there have been a number of studies of the 
one-particle and two-particle behavior of electrons lo- 
calized on quantum dots within a quantum computer 
0, 01 . Starting from the simplest case of two electrons 
in singly-occupied dots in the lowest orbital state, sys- 
tematic generalizations have been introduced and their 
effect on the exchange interaction studied. In particular, 
researchers have analyzed the effect of double occupation 

§0, , higher orbital states 0, 0] , many-electron dots 
, and spin-orbit coupling |^. Here we point out that 
many-body interactions can appear once there are more 
than two electron spins in a computer. In a computer 
with at least three coupled dots containing three elec- 
tron spins, multi-particle exchange processes can change 
the strength of exchange interactions. In a computer 
with four coupled dots containing four electron spins, 
four-body interactions can arise. Such modifications to 
the exchange interactions have not been previously ad- 
dressed in the quantum computing literature, and require 
careful consideration, and possibly removal by error cor- 
rection or avoidance if one intends to scale a 
quantum computer beyond two spins. Their ramifica- 



tions for quantum computation may also extend to (i) 
the encoded universality paradigm, where in efficient im- 
plementations several exchange interactions are turned 
on simultaneously 0; (ii) adiabatic quantum computing 
where the final Hamiltonian for any non-trivial cal- 
culation inevitably includes simultaneous interactions be- 
tween multiple qubits; (iii) fault-tolerant quantum error 
correction, where a higher degree of parallelism trans- 
lates into a lower threshold for fault-tolerant quantum 
computation operations 10] ; (iv) the "one-way" quan- 
tum computer proposal, where all nearest-neighbor in- 
teractions in a cluster of coupled spins are turned on 
simultaneously in order to prepare a many-spins entan- 
gled state [13 ; (v) the search for physical systems with 
intrinsic, topological fault tolerance, where systems with 
four-body interactions have recently been identified as 
having the sought-after properties [l3|. 

Our analysis begins with the derivation of an effec- 
tive Hamiltonian for the electron spins in the quan- 
tum computer. The microscopic Hamiltonian describing 
Coulomb-coupled electrons is 



1=1 i<j 



(1) 



where the confining potential ^(r) contains n energy 
minima, which give rise to the n quantum dots that house 
the electrons. We assume at first that each dot contains 
a single energetically accessible orbital, in accordance 
with the Heitler-London (HL) approximation 14], and 
we label the orbitals with capital letters A,B,... Since 
electrons have two possible spin states, this assumption 
leaves the n-electron system with 2" basis states 

\^{SA, SB,...))^J2 SPPUB ...) \SASB ...)]. 
P 

Here, the first ket \AB . . .) refers to the orbital states of 
the electrons and the second ket \sasb ■ ■ ■) indicates the 
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spin of each electron. The sum runs over all n! permu- 
tations P, where Sp = 1 (—1) if the permutation is even 
(odd), ensuring overall antisymmetry. In every term of 
the sum, the electron in orbital A has spin sa- 

In this 2" dimensional basis, the Hamiltonian ^ takes 
the form of a 2" x 2" Hermitian matrix, specified by 
4" real numbers. One obtains an electron-spin repre- 
sentation of the Hamiltonian matrix by writing it as a 
sum of 4" Hermitian spin matrices of the form (Jij^... = 
(Ji [^A) (g) (Tj (i?) ® • • • each multiplied by a real coefficient 
(there are n factors in Ci j^...). Here, (Jiij)) denotes the 
Pauli matrix tJi acting on the electron in dot p, with 
i = 0,1,2,3 and with ctq equal to the identity matrix. 
This decomposition into spin matrices produces an ef- 
fective electron-spin Hamiltonian that conveniently de- 
scribes the dynamics of n qubits. 

Symmetry considerations can fundamentally constrain 
the form of the electron-spin Hamiltonian. For simplicity, 
we will assume that the quantum dots in the computer 
are arranged in a completely symmetric fashion (i.e., an 
equilateral triangle for 3 electrons, an equilateral tetrahe- 
dron for 4 electrons). We will also neglect spin-orbit cou- 
pling and any external magnetic field. These assumptions 
are introduced for convenience - they are not essential to 
our conclusions, and will be relaxed in a future publica- 
tion _l_5j . They simplify the analysis by implying that the 
effective spin operator Hamiltonian has rotation, inver- 
sion, and exchange symmetry. It can therefore only be a 
function of the magnitude squared — (8^1 + 85 + . . . )^ 
of the total spin 8t- We must have 



distinct eigenvalues of 8|, is [^J -I- 1 (where [§J denotes 
the greatest integer less than ^), so this is the maximum 
number of distinct energy eigenvalues of the Hamiltonian 
The Hamiltonian need only contain this many de- 
grees of freedom to fix all of its eigenvalues, so we can set 
i,„ = for TO > Lf J + 1- We are led to [f J 1 coupled 
linear equations for the non-zero Lm parameters. 

In the three electron case, the total spin can be St = 
1/2 or St = 3/2. We therefore need to solve [|J + 1 = 2 
equations, and it is sufficient to keep only two constants 
Lo and L\ in i/spin, setting Lm>2 = 0. As a convenient 
state with known St = 3/2 we take the normalized state 
1*1^/2) ^ l*(tTT))- Equation ® gives 



E. 
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£3 + 2eo - 3ei 



3/2 - Lo + —Li — — — . 

4 p3 + 2po- 3pi 



(4) 



Here, we have defined 

P3 = {ABC\ABC), 
Pi = {BAC\ABC), 
Po = {CAB\ABC), 



es = {ABC\H\ABC) , 
ei = {BAC\H\ABC) , 
eo = {CAB\ H \ABC) , 



where the subscript i in pi or denotes the num- 
ber of electrons with the same orbital state in bra 
and ket. For the case St = i<t,3 

^(l*(TiT))-|*aTT))) yields 



1/2, using 1*3/^) 



cx 



£■1/2 — 



4 P3-P0' 



(5) 
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Lq + LlSrp + L2(Srp) 



(2) To compute the usual exchange coupling, it is useful to 



where Lq, Li, L2, ■ ■ ■ are real constants with dimensions 
of energy (we take spin to be dimensionless). The con- 
stant Lq is an energy shift. The term proportional to 
Li gives rise to the familiar Heisenberg interaction. Here 
we see that in principle higher order interactions may be 
present in the spin Hamiltonian, starting with a fourth 
order term proportional to ^2- (The constants L„ are 
expected to decrease in magnitude with n but we will 
demonstrate that at least L2 can be physically impor- 
tant in quantum computations.) 

To compute the values of Lo, Li, ^2, • ■ • of the n elec- 
tron effective-spin Hamiltonian, we consider an eigen- 
state l^'s^) of 8y, with known eigenvalue St{St + !)■ 
We compute the expectation value of the effective spin 
Hamiltonian Q in this state, compute the expectation 
value of the spatial Hamiltonian (^), and equate the two: 



*sJ^.pin|*SJ = 



(3) 



This procedure is repeated for all eigenvalues of 8^, thus 
generating a set of linear equations for the parameters 
Lq, ii, ^2, • ■ ■ , in terms of matrix elements of H between 
different orbital states. For n electrons the number of 



rewrite -ffspin as 

L^spin — [Lq + Li ^ 8^) + 2Li ^ 8i • 8j 

A<i<C A<i<j<C 

= if + J(8^.8b + 8^.8c + 8b.8c). (6) 

The constants K ^ Lq + jLi and J — 2Li can be ex- 
pressed in terms of the pi and using Eqs. Q and 
Note that the value of the exchange constant J is deter- 
mined in part by the "three-electron-exchange" terms of 
the form po = {CAB \ABC) and eo = {CAB\ H \ABC). 
Such terms involve a cooperative effect between all three 
electrons and hence cannot be seen in two-electron calcu- 
lations. Thus, the presence of the third electron quanti- 
tatively changes the exchange coupling between the other 
two electrons. 

To compute values for the pi and e^, it is necessary to 
select a specific form for the one-body potential in 
We choose the sample potential 

V(^)-^^^l\^-^\'\^-'B\'\^-C\'\r-D\\ (7) 

This potential has a quadratic minimum at each of 
the vertices of an equilateral tetrahedron A = (0,0,0), 



FIG. 1: Plot of the exchange constant as a function of tunnel- 
ing barrier and Coulomb energy in the case of three mutually 
interacting electrons (dimensionless units). 



FIG. 2: Plot of the change in the exchange constant due to 
3-electron-swap matrix elements. Horizontal axes as in Fig. 
(with their directions flipped). 
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-l,—2l\lj), which become the locations 



(2?y'i, 0,-2/^1), C = {-1^^,1,-21^1), ^^'^ 

of the dots A,B,C,D. The distance between vertices is 
21. We select a potential with four minima so that it 
can be used in the four electron case without modifica- 
tion; the extra minimum does not influence the three- 
electron case in any significant way. At vertex A, we 
define the localized Gaussian state (/)A(r) = (r|^) = 
(^)^^'*exp (-^|r - A|2) which is the ground state 
of the quadratic minimum at that vertex. We define lo- 
calized states similarly for the other vertices. The Gaus- 
sian form makes it possible to obtain all pi and an- 
alytically in terms of the energy fiujo, and the dimen- 
sionless parameters Xb = mijJoP /Ti and Xc = / {nlhuio)- 
Physically, Xf, is the ratio of the tunneling energy barrier 
imtj^P to the harmonic oscillator ground state energy 
Eg = ^hojo, while Xc is the Coulomb energy over Eg. 

In Fig. ^ we plot the exchange-interaction constant J 
in units of huo as a function of Xf, and be. The plot gen- 
erally indicates that J increases as the tunneling barrier 
decreases {xb smaller), an intuitively reasonable result. 
A negative minimum is visible when ~ 5; this may 
reflect the breakdown of the HL approximation in the 
regime of strong interactions. Following Ref. Q, we esti- 
mate that realistic values for GaAs heterostructure single 
dots, Xb ~ 1, ~ 1.5, and hujo ~ 3 meV, do not fall in 
this strong interactions regime. 

Fig. 121 shows the change A J that results from three- 
electron-swap matrix elements eo and po (i.e. A J is the 
change in J = 2Li when eo and po are set to zero in 
expressions (0J and (jS])). Comparing the scales of Figs. 
and 121 one finds that the three-electron swap matrix el- 
ements can have a powerful influence on J. They will 
strongly impact quantum evolution when two-particle 
gates act simultaneously on three qubits (as may arise 
in circumstances (i) - (iii) above). 



In the four electron case, a still more striking effect 
arises. Here, it is possible to have St = 0, St — ^, or 
St — 2, so that one keeps three constants Lq, Li, and 
L2 in -ffspin- It follows immediately that -ffspin includes 
terms of the form L2{Sa'Sb)(Sc 'So) and permutations. 
Unless L2 happens to vanish, the presence of a fourth 
electron introduces a qualitatively new 4-body interaction 
as well as a quantitative change in the exchange coupling 
between the other electrons. 

We evaluate Lq, Li, and L2 using Eq. Extending 
the three-electron deflnitions, we set 



Po = 


{BADC\ABCD), 


eo = 


{BADC\H\ABCD), 


P'o- 


{DABC\ABCD), 


^'0 = 


{DABC\ H \ABCD) , 


Pi = 


{ADBC\ABCD), 


ei = 


{ADBC\H\ABCD), 


P2 = 


{BACD\ABCD), 


£2 = 


{BACD\H\ABCD), 


P4 = 


{ABCD\ABCD), 


£4 = 


{ABCD\H\ABCD). 



A convenient state to use for St = is |5'q) oc 

(i*(nn)> - i^aiit)) - i^arn)) + i^atit))). After 

normalization, this state yields the singlet energy 



En — Ln — 



64 - 4ei + 3eo 
P4 - 4pi + 3po ' 



A convenient state to use for St — I is (x 

(i^-aiti)) + i*(niT)) - i^atti)) - i^-atit))). This 

state, after normalization, yields the triplet energy 



El 



2Li + 4^2 



e4 - 2t2 - ep + 2£o 
P4 - 2p2 - Po + 2po 



Finally, a convenient state to use for S't = 2 is I^Pl) cx 
hl^(TTTt)} • We flnd for the quintet energy 



^2 = io + 6L1 + 36^2 = 



£4 - 6£2 + 8£i -I- 3£o — 6£q 

P4 - 6p2 + 8pi + 3j3o - 6po 
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FIG. 3: Plot of the exchange constant in the case of four 
mutually interacting electrons. Horizontal axes as in Fig.Q 




FIG. 4: Plot of the 4-body exchange constant in the case 
of four mutually interacting electrons. Horizontal axes as in 
Fig. El 

We would like to exhibit interaction constants explic- 
itly in the spin Hamiltonian. Spin-operator identities al- 
low one to convert from the form Q to 

i7,pi„ = + J^S, .8, + J'[(S^-Sb)(Sc-Sz,) 

+ {Sa ■ Sc) {Sb ■ Sd) + (Sa ■ Sd) {Sb ■ Sc)] 

where K = Lq + 3Li -I- ^^2, J = 2Li + 14L2, and 
J' = 8^2- Generically, J' does not vanish, and four-body 
interactions arise. 

We use the potential |7J to estimate the magnitude of 
the effect. The exchange- interaction constant J is plotted 
in Fig. 13 Its form is similar to that of J in the three- 
electron case, with a slightly reduced magnitude. The 
four-body interaction constant, J', appears in Fig. 0] In 
the parameter region of likely experimental interest {xh ^ 
1, Xc ~ 1.5, hujo ~ 3 meV), J'/J ~ —15%, certainly 
large enough to demand attention in computer design. 
Experimentally, four-body terms have been observed in 
•^He 0, and CU4O4 square plaquettes in La2Cu04 0|, 
where J'/J was found to be ^ 27%. 



Finally, to test the accuracy of our HL calculations, 
we have used a Hund-MuUiken (HM) calculation that 
extends the HL basis to include basis states with dou- 
ble occupation: two electrons on a single dot. For in- 
stance, the HM Hilbert space describing three electrons 
on three dots splits into a four-dimensional total-spin 
subspace with St = 3/2 (same as in the HL case), and 
a 16-dimensional St — 1/2 subspace, comprised of two 
degenerate St, 2 = ±1/2 subspaces. For parameters of in- 
terest 01, twelve St — 1/2 states have high energy and 
involve substantial double occupation while the remain- 
ing four 5*7- — 1/2 states form a degenerate ground state 
similar to that of the HF basis. Therefore one recovers 
the HF picture by projecting out the 8 lowest members 
of the HM solution (the four low St = 1/2 and the four 
St — ?>I2 states) to obtain an exchange Hamiltonian 
The corresponding J qualitatively confirms the form of 
Fig. n and Fig. |21 Thus, even when HL results arc not 
necessarily quantitatively precise, substantial many-body 
corrections to the interaction Hamiltonian persist in more 
accurate calculations. In the context of quantum compu- 
tation, these effects could, on one hand, utterly derail an 
algorithm or, on the other, find uses in novel designs. 
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